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ABSTRACT 

In the inner regions of accretion disks around compact objects, the orbital frequency of the gas deviates from 
the local Keplerian value. For long-wavelength modes in this region, the radial epicyclic frequency k is higher 
than the azimuthal frequency 57. This has significant implications for models of the twin kHz QPOs observed 
in many neutron-star sources that traditionally identify the frequencies of the two kHz QPOs with dynamical 
frequencies in the accretion disk. The recognition that the highest frequency in the transition or boundary 
region of the disk is actually the epicyclic frequency also modifies significantly the constraints imposed by the 
observation of high-frequency QPOs on the mass and radius of the compact objects. 
Subject headings: accretion disks — waves 
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1. INTRODUCTION 

The X-ray brightness of an accreting compact object is of- 
ten observed to be modulated quasi-periodically at different 
characteristic frequencies that are comparable to the dynam- 
ical timescale of the central object (see, e.g., van der Klis 
2005). The physical origin of the various types of quasi- 
periodic oscillations (QPOs) in accreting black holes and neu- 
tron stars are still a matter of debate (see, e.g., Psaltis 2004; 
van der Klis 2005). Indeed, different variability models at- 
tribute the observed QPOs to brightness variations generated 
at different dynamical frequencies at a particular radius in the 
accretion disk (e.g., Alpar & Shaham 1985; Miller, Lamb, 
& Psaltis 1998; Stella, Vietri & Morsink 1999; Abramowicz 
et al. 2003) or at the frequencies of wave modes in the disk, 
which are also related to the dynamical frequencies (e.g., Al- 
par et al. 1992; Alpar & Yilmaz 1997; Wagoner 1999; Kato 
2001). 

Despite their differences, most proposed models agree that 
the highest frequency of a large amplitude QPO cannot be 
larger than the frequency of a stable dynamical oscillation in 
the accretion disk. Around a non-rotating black hole or a rel- 
atively compact neutron star of mass M, the highest stable 
dynamical frequency is the azimuthal orbital frequency at the 
radius of the innermost stable circular orbit 
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Requiring this frequency to be larger than the highest ob- 
served QPO frequency imposes an upper bound on the mass 
of the compact object (e.g.. Miller, Lamb, & Psaltis 1998) 
or even provides evidence that the compact object is rapidly 
spinning (e.g., Strohmayer 2001). 

In a realistic accretion disk, the dynamical frequencies of 
oscillations of fluid elements are approximately equal to the 
dynamical frequencies of test particles (such as eq. (Q) only 
away from the boundaries. For example, in the accretion disk 
around a neutron star with a dynamically important magnetic 
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field, the orbital frequency of a fluid element near the so- 
called Alfven radius (see, e.g., Ghosh & Lamb 1991) deviates 
significantly from the local Keplerian frequency because of 
magnetic stresses and viscous stresses (Erkut & Alpar 2005). 
In accretion flows without large scale magnetic fields, such 
as those around a non-magnetic neutron star or a black hole, 
the Maxwell stresses due to turbulent small-scale magnetic 
fields can also affect the dynamical frequencies in the inner 
disk (see, e.g., Hawley & Krolik 2001). Even in the absence 
of any magnetic fields, radiation drag forces can alter signif- 
icantly the dynamical frequencies at the inner regions of ac- 
cretion disks (Miller & Lamb 1996). 

The magnetic and radiation forces alter the dynamical fre- 
quencies predominantly near the inner regions of the Keple- 
rian flows, where the observed QPOs are expected to origi- 
nate. In this Letter, we argue that in realistic accretion disks, 
in which the azimuthal orbital frequency has a maximum at 
some radius outside the surface or horizon of a compact ob- 
ject, the radial epicyclic frequency at a comparable radius is 
the highest dynamical frequency in the system. Our result has 
significant implications for models of QPOs in accreting com- 
pact objects and for the constraints imposed on the masses 
and spins of the compact objects by the observation of high- 
frequency QPOs. 

2. THE HIGHEST DYNAMICAL FREQUENCY 

In order to discuss the relative order of different dynami- 
cal frequencies in an accretion disk independently of the de- 
tails of the additional forces that affect the motion of fluid 
elements, we use the simple but transparent derivation of the 
equations of motion discussed by Hill (1878). These equa- 
tions are valid for the motion of test particles in a frame ro- 
tating with an angular velocity flo, in the presence of addi- 
tional forces. This set of equations has been used by Balbus & 
Hawley 1992 to elucidate the physics behind the magnetoro- 
tational instability in MHD accretion disks (see also Pessah & 
Psaltis 2005 for a discussion of a more general case of MHD 
instabilities). 

We start with particles in a stable circular orbit at a radius ro 
and with an angular velocity ilo = il(ro), i.e., r = ro, (f) = Hot, 
where n(r) is not necessarily the Keplerian angular velocity. 
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Introducing small perturbations around the circular orbit, i.e., 
r=ro+x and = fio? +y/ro, to first order in x and y, we obtain 
the Hill equations 



x—2flQy = -x 
y+2r2oi = /v 



(2) 
(3) 



where and /, represent the x- and y- components of the 
sum of the perturbations of all forces besides gravity that are 
acting on the particles. For small displacements away from 
the stable equilibrium orbit, f^ and /,, are linear in x and y 
with negative first derivatives, and hence 
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Looking for solutions of the form > 
persion relation 
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we obtain the dis- 



Of. 



dx 



dn^ 

dlnr 







dfy 



dy 



dx 







Lo —liilo 



dfy 



dfy 



dy 



where 



2 _ 



dx 



dn^ 



dx 
df. 



df. 



dy 



dy 



^0, 



d\nr 
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When non-gravitational forces are absent, the solutions are 
stable modes with lo = Q and lo = n. As can be verified eas- 
ily, a; = K is the frequency associated with the mode of ra- 
dial oscillations. As azimuthal perturbations simply shift the 
phase angle cj) among equivalent phases on the stable circular 
orbit ("spontaneous symmetry breaking"), there is no restor- 
ing force, and the solution corresponding to azimuthal oscilla- 
tions is simply a; = 0. Indeed, after an azimuthal perturbation 
the particle would proceed from its shifted phase on the same 
circular orbit, at the original orbital frequency fio- This is an 
instance of the Goldstone theorem that zero frequency modes 
are associated with spontaneous symmetry breaking. The 
physical reason for the nonzero restoring force and nonzero 
frequency lo = n for radial perturbations is rotation. The long 
range correlations introduced by the presence of rotation, i.e., 
the Coriolis force and the gradient of the centrifugal potential 
provide the restoring force in the Hill equations 0. Modes ef- 
fected by long range forces always have nonzero frequencies, 
as exemplified by the plasma frequency due to the Coulomb 
force, the Jeans mass associated with gravitational instabili- 
ties, and the finite mass of Higgs bosons. 

The effect of short-range forces can be traced in the realistic 
limit 
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The LO ~ K modes are now quasi-stable with TI{lo) ~ k and 
T{lo) ~ 0{\df_x.y/ dx]^/ k) <C TI{lo). There is still an exact lo = 
mode, due to the axial symmetry and a decaying mode with 



TZ{lo) = and T(lo) ^ —0(k). This decay corresponds to the 
effect of the shear on the initial perturbations. 

When the disk is treated as a continuous medium and its hy- 
drodynamic or magnetohydrodynamic wave modes are stud- 
ied, the analogues of restoring forces on perturbations of sin- 
gle particle orbits in the Hill equations are spatial derivatives 
of the forces on a fluid element. These restoring force terms 

contain the product (Va ■ k)^ or (csk)^ where k is the wave vec- 
tor, and Va and are the Alfven velocity and sound speed 
respectively (see, e.g., Alpar et al. 1992; Balbus & Hawley 
1992; Alpar & Yilmaz 1997; Kato 2001; Pessah & Psaltis 
2005). Conditions analogous to (|8} describe global modes of 
long wavelength. These modes are expected to be the ones 
that produce the largest amplitude of brightness modulation. 
In the limit of small wavenumber k, the frequencies of os- 
cillations have values close to k and ft. For modes with m 
azimuthal nodes, frequencies comparable to k - mft will ap- 
pear; the dominant global modes will be those with m = and 
m= 1, which represents the beat between the radial epicyclic 
frequency and the disk rotation. This fundamental effect of 
rotation can be seen in all dispersion relations of waves in ro- 
tating fluids and accretion disks, in models of varying degrees 
of complexity (see, e.g., Chandrasekhar 1961; Papaloizou & 
Pringle 1984). We emphasize here that not all oscillations de- 
scribed by the Hill equations are stable. This is of particular 
relevance to MHD disks that are subject to the magnetorota- 
tional instability, especially at the limit of small wavenum- 
bers. 

Where magnetic or radiation forces significantly affect the 
orbital motion of gas elements, the azimuthal frequency in the 
disk differs from the local Keplerian value. As an example, 
for accretion onto a rotating and magnetic neutron star, the 
local azimuthal frequency il{r) is lower than the Keplerian 
frequency inside some characteristic radius r2 comparable to 
the Alfven radius. It has a maximum ilmax at r < r2 and then 
decreases to match the rotation frequency of the star 51* at 
a radius ri , which marks the disk-magnetosphere boundary. 
The simplest mathematical model employs a quadratic form 
for n{r). 
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where M is the mass of the star. Matching these two expres- 
sions, as well as their derivatives, at r2 relates the maximum 
frequency and its location to the transition region parameters 
ri, r2\ and the neutron-star rotation rate fi*: 

_ r>,/r>2-(9/16)(ri/r2-7/3)^ ^ _1 ( ^ 3 r ^^ ^ 

f}.M + (3/2)(n/r2-5/3) '''-4lv Vr^j^'^' 

(10) 

This model quite accurately represents the Q,{r) curves that 
meet a boundary condition at the stellar magnetosphere and 
asymptotically join the Keplerian curve within a transition re- 
gion in models of accreting neutron stars with dynamically 
important magnetic fields Erkut & Alpar 2005. The model 
of the azimuthal frequency profile described by equation (|9j 
is shown in Figure 1, for ^2/^* = 1-1, ^max/^* = 1-2, and 
fi/f'2= 0.84, together with the radial epicyclic frequency cal- 
culated according to equation (|6|l. 

At radii r < r2, where the azimuthal frequency is lower 
than the Keplerian frequency, the radial epicyclic frequency 
is the larger of the two dynamical frequencies. At the radius 
r, where the azimuthal frequency has a maximum, k = 2ilniax 
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Fig. 1 . — The radial profile of the azimuthal (f2) and radial epicyclic (k) 
frequencies of the simple model discussed in the text, for = 11 and 

ri/r2 = 0.4. The dashed line shows the difference k — Q. At radii r > ri 
the azimuthal and radial epicyclic frequencies are equal to the Keplerian fre- 
quency. 



and the ratio n/fl increases to even larger values in the region 
ri < r <r, where the azimuthal frequency is decreasing with 
decreasing radius. For a wide range of radii inside r2, the 
radial epicyclic frequency is larger than the local Keplerian 
frequency. As a result, if low-^ modes in this transition region 
are responsible for the observed QPOs in accreting neutron 
stars and black holes, as envisioned by most models, then an 
upper bound of the QPO frequencies may be set not by the 
maximum azimuthal frequency but by the maximum radial 
epicyclic frequency in the region. 

3. DISCUSSION 

In §2, we argued that, in the inner regions of accretion disks 
around compact objects, magnetic, radiation, pressure, and 
viscous forces may become comparable to the gravitational 
force, alter the orbital frequencies of fluid elements and lead 
to radial epicyclic frequencies in excess of the orbital frequen- 
cies. Here, we discuss the implications of our results for mod- 
els of the kHz QPOs observed from many accreting neutron- 
star sources (for a review see van der Klis 2005). These are 
pairs of QPOs with frequencies comparable to a kHz that vary 
on timescales longer than any of the dynamical timescales in 
the vicinity of the neutron stars. 

As the frequencies of the kHz QPOs vary, they follow a 
number of intriguing patterns. In all sources for which the 
spin frequency of the neutron star is known, either via ob- 
servations of X-ray pulses or burst oscillations, the differ- 
ence frequency between the two kHz QPOs was shown to be 
comparable to the neutron star spin frequency or to half its 
value. This property gave rise to the beat-frequency models 
of kHz QPOs (see Strohmayer et al. 1996; Miller et al. 1998; 
Chakrabarty et al. 2003). In sources for which the frequencies 
of the two kHz QPOs were measured with high accuracy, such 
as Sco X-1, it was shown that the two frequencies followed a 
quadratic correlation (Psaltis et al. 1998; Psaltis, Belloni, & 
van der Klis 1999). Together with the correlation between 
the frequencies of the kHz QPOs and other low-frequency 
QPOs observed simultaneously in the same sources, this pro- 
vided support for the relativistic precession models (Stella & 
Vietri 1999; Stefla, Vietri, & Morsink 1999; Psaltis & Nor- 
man 2000). Finally, the ratios of the frequencies of the two 
kHz QPOs are roughly comparable to the ratio 2/3 and this 



gave rise to models that rely on resonances between modes 
(Abramowicz et al. 2003). It is important to note here that, 
mathematically speaking, not all three of the above patterns 
can characterize simultaneously the frequency correlations of 
the kHz QPOs and, in fact, none of them is valid to within 
the measurement errors of the QPO frequencies. However, all 
three of them can be shown to be able to describe the data 
roughly, typically to within 30% for the first and last alterna- 
tive and to within 5% for the quadratic correlation (see Psaltis 
et al. 1998 for a detailed discussion of the statistical signifi- 
cance of various correlations). 

In all the models of the kHz QPOs mentioned above, one 
or both of the observed QPO frequencies are interpreted as 
dynamical frequencies (azimuthal, radial, or vertical) at dif- 
ferent characteristic radii in the accretion flows. For each of 
these models, these dynamical frequencies have been calcu- 
lated so far assuming that the only force that affects the mo- 
tion of fluid elements is gravity, even though there is always 
an implicit assumption that some additional physical mecha- 
nism picks the characteristic radius where these QPOs origi- 
nate. As we have shown, the relative ordering of the frequen- 
cies for long-wavelength modes depends on the radial pro- 
file of the azimuthal frequency. Therefore, the identification 
of observed QPOs with particular dynamical frequencies may 
not be self consistent in any of the above models. 

For most realistic models of the azimuthal frequency, e.g., 
for the quadratic model discussed in §2 (eq. (9l), the largest 
of the dynamical frequencies is k and the smallest is 51, in the 
inner part of the transition region where dfl/dr > (see also 
Fig. 1). Wavepackets in the disk can modulate the accretion 
flow onto the compact object at the radial epicyclic frequency 
band ^ k, at the orbital frequency band ^ ft, or at the beat 
K-ft of the radial and orbital motions. It is clear that k > 2i7 
and therefore k - 57 > O in the region dft /dr > 0, where the 
order of the frequencies is k > K-n> ft. 

How does one identify the observed upper and lower kHz 
QPO frequencies, hereafter 1^2 and vi, respectively, with two 
of these three possible frequencies? There are two simple 
trends that all kHz QPOs are observed to obey: (i) the two 
kHz QPO frequencies always increase or decrease together, 
and (ii) their difference Ai/ = V2-v\ decreases as both fre- 
quencies increase. The only pair of frequencies that satisfies 
these two conditions are k and K-f2. Indeed, if we associate 
these two frequencies with the upper and lower kilohertz QPO 
frequencies, i.e., k, = Ittvi and n-Vl = 2ttv\, they will satisfy 
the trend (ii) if dil/dK < 0, which is equivalent to 

diogndn/dr) 

diogr 

This condition holds as long as the transition region size is 
about a fifth of the disk inner radius or less. 

Perhaps the most interesting application of kilohertz QPO 
models is the constraints they provide on the mass-radius rela- 
tion of neutron stars (Miller et al. 1998; see also discussion af- 
ter Eq. [1]). The interpretation of the highest QPO frequency 
as the epicyclic rather than the Keplerian azimuthal frequency 
modifies the bounds obtained from the upper kilohertz fre- 
quency. To estimate the effect of the phenomena discussed 
here on the bounds in the mass-radius plane, we use the simple 
model (|9jl. In this model, the maximum epicyclic frequency, 
to a first approximation, is equal to the radial epicyclic fre- 
quency at radius r and hence 

27riy2,max ^ K(f) = 2^^.^^ . (12) 
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Fig. 2. — Constraints on the mass-radius plane of neutron stars imposed 
by the observation of the highest frequency QPO so far (1330 Hz, se e va n 
Straaten et al. 2000) for different values of the parameter / (see, eq. 1141 1. 
The solid lines show the predictions of representative equations of state for 
neutron-star matter (see Lattimer & Prakash 2001). 



The constraints the kHz QPOs impose on the mass and 
the radius of a neutron star are obtained from two require- 
ments. First, the neutron star radius tns must be less than a 
representative inner disk radius that is associated with upper 
kilohertz QPO frequency, which was interpreted as a Keple- 
rian frequency in earher appHcations (see van der Klis 2005). 
Recognition of the epicyclic frequency as the highest fre- 
quency modifies this constraint to 
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depends on the width of the transition region. The second 
constraint simply states that the radius of the last stable orbit 
is less than the disk radius associated with the upper kilohertz 
QPO, risco = dGM/c^ < r. This leads to 
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if the upper kilohertz QPO is the radial epicyclic frequency 
rather than the Keplerian frequency. For an infinitely nar- 
row boundary layer, f = 2, while for r/r2 = 0.9 and 0.8, f = 
1.84 and 1.65 respectively. The resulting modified constraints 
are shown in Figure 2. We conclude that a careful interpreta- 
tion of the kHz QPO frequencies makes the constraints on the 
masses and radii of neutron stars much less stringent. 

We note that in general relativity, the radial epicyclic fre- 
quency for test particle orbits deviates from the Keplerian 
form at radii close to risco- Unlike the Newtonian case, in 
general relativity, radial and vertical epicyclic frequencies do 
not coincide with the azimuthal frequency. Simple derivations 
of epicyclic oscillation frequencies in the (Newtonian, as well 
as) general relativistic cases are provided by Abramowicz & 
Kluzniak (2004), who give expressions for the dynamical fre- 
quencies in terms of the metric. General relativistic expres- 
sions for the radial epicyclic frequency incorporating fluid ef- 
fects and appropriate formulations of the metric around the 
neutron star will be the subject of future work. By continu- 
ity from the Newtonian treatment here, we anticipate that the 
highest frequency will turn out to be the radial epicyclic fre- 
quency at radii ri,r > risco by factors of a few. 
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